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1 Introduction 



A comprehensive treatment of the theory of Grobner bases is far beyond the 
scope of this tutorial. The reader is reffered to Q], [2J, [5] for a more detailed 
discussion. Grobner bases offers a powerful algorithmic criterion for the exis- 
tence of solutions to a system of algebraic equations over algebraically closed 
fields. While the theory is rich and builds on the theory of Ideals and Varieties, 
the discussion in the persent tutorial will remain as elementary as possible. The 
tutorial focuses on implementation aspects of Grobner bases computation via 
the Buchberger algorithm. The Sage |S + llj specific implementation discussed 
here is not meant to rival the efficient implementations available on Singular, 
Maxima and other Computer Algebra Systems(CAS). The tutorial merely at- 
tempts to provide an overview of the implementation details of the Buchberger 
algorithm. It must be pointed out that there are many other algorithms besides 
the Buchberger algorithm for computing Grobner bases, however for simplicity 
we restrict our attention to the Buchberger Algorithm. The reader can contact 
the author Q , to obtain the I^TjrjX source or to suggest possible corrections 
The problem that the tutorial discusses is the determination of existence of 
solutions to a system of polynomial equations of the form. 

fi(xi, ■■■ ,x n ) = ■•■ = f m (xi, ■■■ ,x n ) = 0. (1) 

where 

VI < k < m fk(xi,- ■ ■ ,x n ) e C[xi, ■ ■ ■ ,x n ]. 

The determination of existence of solution can be established algorithmically 
via Grobner bases computation. 



2 Initial Set up 

For simplicity we will be considering systems of polynomial equations with at 
most 10 variables. 
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# Defining the variables 

varCxO, xl, x2, x3, x4, x5, x6, x7, x8, x9') 

# The prime assignment 



p = 


= Primes () 






pO 


= P .unrank(O) ; 


Pi 


= P.unrank(l) 


P2 


= P.unrank(2) ; 


p3 


= P.unrank(3) 


p4 


= P.unrank(4) ; 


p5 


= P.unrank(5) 


p6 


= P . unrank(6) ; 


P7 


= P.unrank(7) 


p8 


= P.unrank(8) ; 


p9 


= P.unrankO) 



3 Multivariate leading term 

The multivariate leading term functions determines the leading term of an in- 
put polynomial. The leading term for a given polynomial is identified using a 
monomial ordering. We use in the tutorial a monomial ordering proposed by 
the author which follows from the fundamental theorem of arithmetics. The 
proposed ordering has the advantage of considerably simplifying aspects of the 
implementation. 

def multivariate_leading_term(f ) : 

# Expression used for specifying the 

# type of the operation, 
add = xO+xl 

mul = xO*xl 
xpo = xO~xl 
est = 2 

if (f . operator () == add. operator ()) : 

# Collecting the terms 
L = f.operandsO 

# Collecting the terms striped from their 

# coefficients 
L_strpd = listO 

for i in range (len(L) ) : 

if (L [i] ). arguments () != () : 

if ( (L [i] ). operator () == mul . operator () or \ 
(L [i] ). operator () == xpo . operator ()) : 
cst_fctr = 1 

lst_i = (L [i] ). operands () 
for j in range (len(lst_i) ) : 

if (lst_i [j] ). arguments () == () : 
cst_fctr = lst_i[j] 
L_strpd. append ( (expand(L [i] /cst_f ctr) , i) ) 
elif (L[i] == xO or L[i] == xl or L[i] == x2 or \ 
L[i] == x3 or L[i] == x4 or L[i] == x5 or L[i] == x6 or \ 
L[i] == x7 or L[i] == x8 or L[i] == x9 ): 
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L_strpd. append ((L[i] ,i)) 
# Storing the integer and the index associated with the term 
tmp_value = (L_strpd [0] [0] ). substitute (x0=p0 ,xl=pl ,x2=p2 , \ 
x3=p3 , x4=p4 , x5=p5 , x6=p6 , x7=p7 , x8=p8 , x9=p9) 
idx = L_strpd[0] [1] 
for k in range (len(L_strpd) ) : 

if (L_strpd[k] [0] ) .substitute (x0=p0,xl=pl,x2=p2,x3=p3,\ 
X4=p4,x5=p5,x6=p6,x7=p7,x8=p8,x9=p9) > tmp_value : 

tmp_value = (L_strpd [k] [0] ). substitute (x0=p0 ,xl=pl , \ 
x2=p2 , x3=p3 , x4=p4 , x5=p5 , x6=p6 , x7=p7 , x8=p8 , x9=p9) 
idx = L_strpd[k] [1] 
return L [idx] 
elif (f . operator () == mul . operator ()) : 
return f 

elif (f . operator () == xpo . operator ()) : 

return f 
else : 

return f 



4 The multivariate division function 

The multivariate division function is at the heart of the Buchberger algorithm 
for computing Grobner bases. This routine generalizes the familiar long division 
algorithm for single variable polynomials two non trivial ways. First the algo- 
rithm extends the single variable long division algorithm to polynomials with 
multiple variables. Second the algorithm allows for the division of single multi- 
variate polynomial by a finite list of multivariate polynomials. It goes without 
saying that the implementation relies quite heavily on the monomial ordering . 
We provide here both the pseudocode description of the algorithm as found in 
[T] and our proposed Sage implementation. 
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1. Input: {fk}!< k < s , / 

2. Output: {afc}i<fc< s , r 

3. Initialization^ <— 0} 1<fc<s ; r ; p 4— f 

4. WHILE p/ODO 

(a) i <- 1 

(b) divisionoccured <— false 

(c) WHILE i < s divisionoccured =false DO 

i. IF LT(fi) divides LT(p) THEN 
A „.<_„. + L7 » 



C. divisionoccured -s— true 
end 
ii. ELSE 

A. + l 
end 

end 

(d) IF divisionoccured = FALSE 

i. r <— r + LT(p) 

ii. p <— p — LT(jp) 

end 

end 
5. End 
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The implementation of the algorithm described above is also particularly 
insightfull. 

def multivariate_division(f , List): 

# Initializing the output List 
L = listO 

for j in range(len(List)+l) : 
L . append (0) 

# Initializing the Polynomial 
P = f 

while ( p != ) : 
i = 

division_occured = 

while (i in range (len(List) ) ) and (division_occured == 0): 
if List[i] == 0: 

i = i+1 
else : 

# Getting the Leading term of fi 

Lt_fi = multivariate_leading_term(List [i] ) 
#Getting the leading Monomial of fi 

Lm_fi = Lt_fi/Lt_fi. substitute (x0=l,xl=l,x2=l,x3=l,\ 
x4=l ,x5=l ,x6=l ,x7=l ,x8=l ,x9=l) 

# Getting the Leading term of p 
Lt_p = multivariate_leading_term(p) 
#Getting the leading Monomial of p 

Lm_p = Lt_p/Lt_p. substitute (x0=l,xl=l,x2=l,x3=l,\ 
x4=l ,x5=l ,x6=l ,x7=l ,x8=l ,x9=l) 

m_p = Lm_p . substitute (x0=p0,xl=pl ,x2=p2 ,x3=p3 , \ 
x4=p4 , x5=p5 , x6=p6 , x7=p7 , x8=p8 , x9=p9) 

m_fi = Lm_fi . substitute (x0=p0 ,xl=pl ,x2=p2 ,x3=p3, \ 
x4=p4 , x5=p5 , x6=p6 , x7=p7 , x8=p8 , x9=p9) 

if (gcd(m_p, m_fi) == m_fi or gcd(m_p, m_fi) == -m_f i) : 
L[i] = expand(L[i] + Lt_p/Lt_fi) 
p = expand (p - List [i] * (Lt_p/Lt_f i) ) 
division_occured = 1 
else : 

i = i+1 
if (division_occured == 0) : 

L[len(List)] = L[len(List)] + multivariate_leading_term(p) 
p = p - multivariate_leading_term(p) 

return L 
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5 The multivariate least common multiple 

The multivariate least common multiple function determines the monomial least 
common multiple(LCM) for an input pair of monomials. Our proposed mono- 
mial ordering reduces the monomial LCM computation to the familiar integer 
LCM computation. 

def multivariate_monomial_lcm(tl , t2) : 

# These 2 lines of code get rid of the coefficient of the leading terms 
ml = tl/tl . substitute (xO=l ,xl=l ,x2=l ,x3=l ,x4=l ,x5=l ,x6=l ,x7=l ,x8=l ,x9=l) 
m2 = t2/t2 . substitute (xO=l ,xl=l ,x2=l ,x3=l ,x4=l ,x5=l ,x6=l ,x7=l ,x8=l ,x9=l) 

# The following computes the 1cm value associated with the monomial we seek 
monomial_lcm_value = 1cm (ml . substitute (xO=pO ,xl=pl ,x2=p2 ,x3=p3 , \ 

x4=p4 ,x5=p5 ,x6=p6 ,x7=p7 ,x8=p8 ,x9=p9) , m2 . substitute (xO=pO ,xl=pl ,x2=p2 , \ 
x3=p3 , x4=p4 , x5=p5 , x6=p6 , x7=p7 , x8=p8 , x9=p9) ) 



# The next section of line of codes recovers the 

# monomial in question from the computed 1cm integer. 
prime_f actors = factor (monomial_lcm_value) 
factor_list = list (prime_f actors) 

m = 1 

for i in range (len(factor_list) ) : 
tmp_list = list (f actor_list [i] ) 
if (tmp_list [0]==p0) : 

m = m*xO~Integer (tmp_list [1] ) 
elif (tmp_list [0]==pl) : 

m = m*xl~Integer (tmp_list [1] ) 
elif (tmp_list [0]==p2) : 

m = m*x2~Integer (tmp_list [1] ) 
elif (tmp_list [0]==p3) : 

m = m*x3~Integer (tmp_list [1] ) 
elif (tmp_list [0] ==p4) : 

m = m*x4~Integer (tmp_list [1] ) 
elif (tmp_list [0] ==p5) : 

m = m*x5~Integer (tmp_list [1] ) 
elif (tmp_list [0] ==p6) : 

m = m*x6~Integer (tmp_list [1] ) 
elif (tmp_list [0] ==p7) : 

m = m*x7~Integer (tmp_list [1] ) 
elif (tmp_list [0]==p8) : 

m = m*x8~Integer (tmp_list [1] ) 
elif (tmp_list [0] ==p9) : 

m = m*x9~Integer (tmp_list [1] ) 

return m 
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6 The multivariate substraction polynomial 



The multivariate subtraction polynomial function is an important subroutine 
for the Buchberger Algorithm. Its inner working are quite reminiscent of the 
elimination procedure in Gaussian elimination and it's implementation is just 
as straight forward. 

def multivariate_S_polynomials (List) : 
L = listO 

for i in range (len(List) -1) : 

for j in range (i+1 , len(List) ) : 

Lt_fi = multivariate_leading_term(List [i] ) 
Lt_fj = multivariate_leading_term(List [j] ) 
monomial_lcm = multivariate_monomial_lcm(Lt_f i , Lt_f j) 
Sij = expand(List [i] * (monomial_lcm/Lt_f i) - \ 
List [j] * (monomial_lcm/Lt_f j ) ) 
L . append(Si j ) 

return L 



7 multivariate reduction 

The multivariate reduction function controls the bounds on the degree of the 
polynomials to be adjoined to the Grobner bases. The reduction is achieved by 
dividing the polynomial by the greatest common divisor(GCD) of all the terms 
making up the polynomial. It also follows from our proposed ordering that the 
monomial GCD also reduces to the familiar integer GCD computation. 
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def multivariate_reduce_polynomial (f ) : 

# Checks to see if there is a constant term 

# in which case no reduction is needed 
if f . substitute (x0=0 , xl=0 , x2=0 , x3=0 , \ 

x4=0,x5=0,x6=0,x7=0,x8=0,x9=0) != : 

return f 
elif f == 0: 

return f 
else : 

L = list (f . iterator () ) 



# The next piece of code determines i 

# if the expression is a monomial 
prd = 1 

for j in range (len(L) ) : 

prd = prd*L[j] 
if (prd == f ) : 

return 1 

elif ((len(L) == 2) and (L[0]~L[1] == f)): 
return f 



for i in range (len(L) ) : 

# The next line of code gets rid of 

# the coefficients in the list 

L[i] = L[i]/(L[i] .substitute(x0=l,xl=l,x2=l,x3=l,\ 
x4=l ,x5=l ,x6=l ,x7=l ,x8=l ,x9=l)) 

L[i] = Integer (L [i] . substitute (x0=p0 ,xl=pl ,x2=p2 , \ 
x3=p3 , x4=p4 , x5=p5 , x6=p6 , x7=p7 , x8=p8 , x9=p9) ) 



# Computing the greatest common divisior 
cmn_fctr = gcd(L) 



if cmn_fctr == 1 : 
return f 
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else : 

# The next section of line of codes 

# in question from the computed gcd 
prime_f actors = factor (cmn_fctr) 
factor_list = list (prime_f actors) 

m = 1 

for i in range (len(factor_list) ) : 
tmp_list = list (f actor_list [i] ) 
if (tmp_list[0]==pO) : 

m = m*xO~ Integer (tmp_list [1 
elif (tmp_list [0]==pl) : 

m = m*xl~ Integer (tmp_list [1 
elif (tmp_list [0]==p2) : 

m = m*x2~Integer (tmp_list [1 
elif (tmp_list [0] ==p3) : 

m = m*x3~Integer (tmp_list [1 
elif (tmp_list [0] ==p4) : 

m = m*x4~Integer (tmp_list [1 
elif (tmp_list [0]==p5) : 

m = m*x5~Integer (tmp_list [1 
elif (tmp_list [0]==p6) : 

m = m*x6~Integer (tmp_list [1 
elif (tmp_list [0] ==p7) : 

m = m*x7~Integer (tmp_list [1 
elif (tmp_list [0] ==p8) : 

m = m*x8~Integer (tmp_list [1 
elif (tmp_list [0]==p9) : 

m = m*x9~Integer (tmp_list [1 
g = expand (f/m) 
return g 



recover the monomial 
integer . 



) 
) 
) 
) 
) 
) 
) 
) 
) 
) 
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8 Implementation of the Buchberger Algorithm 



In summary the key components required for the implementation of the Buch- 
berger Algorithm are: 

- The multivariate leading term 

- The multivariate division 

- The multivariate least common multiple 

- The multivariate Substraction polynomial. 

Let us put the pieces together to implement the Buchberger algorithm for com- 
puting Grdbner bases. 

A Grobner bases G of an ideal 1 generated by the set of polynomials {/fc}i<fc<m C 
C[xxj ' ' ' i x io] over the field C is characterised by any one of the following prop- 
erties, stated relative to some monomial order: 

- The ideal given by the leading terms of polynomials in I is itself generated by 
the leading terms of the basis G; 

- The leading term of any polynomial in I is divisible by the leading term of 
some polynomial in the basis G; 

- multivariate division of any polynomial in the polynomial ring C[xi, • • • , Xio] 
by G gives a unique remainder; 

- multivariate division of any polynomial in the ideal I by G gives 0. 

We provide here both the pseudocode description of the algorithm as found 
in pQ and our proposed Sage implementation. 

The Main Theorem and the Buchberger Algorithm 

Let {/fc}i<fc<j ^ be a polynomial ideal. Then a Groebner basis for I can be 
constructed in finitely steps. 

Here is the Buchberger algorithm: 

1. Input: F = fx,-- - ,f s 

2. Output: a Groebner basis G = (g± ■ ■ ■ ,gt) for 1, with F C G 

3. G^F 

4. repeat 

(a) G' <- G 

(b) for each pair {p,q},p ^ q G G do 



i. S^S(p,q) 

ii. if S ^ then G <- G U {S} 
end 

until G=G' 
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Q> 

Where S (p, q) denotes the remainder on division of S (p, q) by the ordered set 
of elements of G' . 

8.0.1 Reduced Grobner bases 

A reduced Grobner bases for a polynomial Ideal I is a Grobner bases for I such 
that : 

1- adeg(f) = 1 for all p E G 

2. For all p 6 G no monomials of p lies in {(LT (G — p))} 
where LT denote the leading term. 

8.0.2 Sage implementation of the Algorithm 

The initialization part which specifies the following system of algebraic equations 
/i := 2x x 2 + 4x1X2 - 6 = / 2 := x\ - x 2 = / 3 := x\-xi = / 4 := Xq - x Q = 

# The Ideal of interest is generated 

# by the following polynomials 

fl = expand((x0+2*xl)*(2*x2)) - 6 
f2 = x2~2-x2 
f3 = xl~2-xl 
f4 = x0~2-x0 

# Generators of the Ideal 

# for solving f l=f 2=f 3=f 4=0 
I = listO 

I . append (fl) 
I. append (f 2) 
I . append(f 3) 
I . append (f 4) 

# Initialization step 
I_curr = listO 

for i in range (len(I) ) : 
I_curr . appendd [i] ) 

l_old = 

l_new = len(I_curr) 
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The main loop for the Algorithm corresponds to the following 

# Boolean variable tracking contradictions 
finished = 

while l_old != l_new and finished == 0: 

# Computes the single pass of the substraction polynomials 
S = multivariate_S_polynomials (I_curr) 
print '\n\n The subtraction polynomials yield' 
for i in range (len(S) ) : 

S[i] = multivariate_reduce_polynomial (S [i] ) 
print >S[>,i,']= ' ,S[i] 

# The instruction bellow is the lazy way of getting rid of the duplicates. 
St = Set(S) 
S = list (St) 

#Recording the size of the Ideal generator set before the division 
l_old = len(I_curr) 
for i in range (len(S) ) : 

tmp_list = multivariate_division(S [i] , I_curr) 
if tmp_list [len(tmp_list) -1] ! =0 : 

I_curr . append (tmp_list [len(tmp_list) -1] ) 

# Printing the result of the first pass of the Buchberger algorithm, 
print '\n\n The Current generator for the Ideal is given by' 
for i in range (len(I_curr) ) : 

print I_curr[i] 

if I_curr[i] == I_curr [i] . substitute (x0=0,xl=0,x2=0,x3=0 , \ 
0,x5=0,x6=0,x7=0,x8=0,x9=0) and I_curr[i] . substitute (x0=0 ,xl=0 , \ 
0,x3=0,x4=0,x5=0,x6=0,x7=0,x8=0,x9=0) != : 
finished = 1 

#recording the size of the generator set after the division 
l_new = len(I_curr) 

I = I_curr 
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9 Running the Script 

For illustration we compute Grobner bases for the Ideal induced by the system 
of equations 

fi := 2x0X2+4x1X2-6 = f 2 := x%-x 2 = /3 := x\-x\ = f 4 := xl-x = 

The output of the computation is the list of size 7 which corresponds to the 
sought after Grobner bases G. 

Gq = 2 x x 2 + 4 xix 2 — 6 



Gi = 




■''2 


G 2 = 


xl- 


■I'l 


G 3 = 


2 

Xq — 


■i'O 


G 4 = 


-x H 


" Xl 


G 5 = 


-x H 


- x 2 


Gq = - 


3 


3 

+ 2 



The criteria for existence of solution to a system of polynomial equations is 
therefore summarized as follows: 

if the inputed polynomial admitted no common root then the Grobner bases 
would include the constant polynomial which constitutes a certificate for non 
existence of solution to the given system of algebraic equations. 

fl(xi, ■■■ ,X n ) = ■■• = f m (xi, ■■■ ,X n ) = 0. 
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